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The structure and the statistical fluctuations of interfaces between coexisting phases in the 
Asakura-Oosawa (AO) model for a colloid-polymer mixture are analyzed by extensive Monte Carlo 
simulations. We make use of a recently developed grand canonical cluster move with an additional 
constraint stabilizing the existence of two interfaces in the (rectangular) box that is simulated. 
Choosing very large systems, of size L x L x D with L = 60 and D = 120, measured in units of 
the colloid radius, the spectrum of capillary wave-type interfacial excitations is analyzed in detail. 
The local position of the interface is defined in terms of a (local) Gibbs surface concept. For small 
wavevectors capillary wave theory is verified quantitatively, while for larger wavevectors pronounced 
deviations show up. For wavevectors that correspond to the typical distance between colloids in the 
colloid-rich phase, the interfacial fluctuations exhibit the same structure as observed in the bulk 
structure factor. When one analyzes the data in terms of the concept of a wavevector-dependent 
interfacial tension, a monotonous decrease of this quantity with increasing wavevector is found. 
Limitations of our analysis are critically discussed. 

PACS numbers: 61.20. Ja,64.75.+g 



I. INTRODUCTION 

By adding non-adsorbing polymers to a colloidal 
suspension, phase separation may be induced. This 
leads to the formation of two coexisting phases, one 
with high colloid density and one with low colloid 
density, separated by an interface. The interface 
is not flat, but subject to thermally driven den- 
sity fluctuations known as capillary waves. Cap- 
illary waves were first predicted by Smoluchowsk: 1 , 
and have since then been studied using light and X- 
ray diffractio n 2 ! 3 ! 4 ! 5 , theoretical method o 6 ! 7 ! 8 ! 9 ! 10 ' 11 , and 
computer simulationsi^iiiii 4 * 1 ^!^!. Recently, capillary 
waves were even visualized directly, close to single parti- 
cle resolution, in a colloid-polymer mixture^. 

Unlike critical fluctuations, capillary waves survive 
deep into the two-phase region of the phase diagram, and 
strongly influence all quantities that depend on transver- 
sal degrees of freedom. For example, the capillary contri- 
bution to the width of the interface diverges logarithmi- 
cally with the lateral system size. As a consequence, the 
apparent width of the interface depends on the length 
scale on which it is studied. This effect is well described 
by capillary wave theory (CWT)^, which consequently 
plays an important role in the modern treatment of in- 
terfaces. The key assumption is that a smooth local in- 
terface position can be defined, see Fig. ^ Clearly, away 
from any critical point, and for long wavelength capillary 
variations, this assumption is reasonable. However, at 
short wavelengths, say of the order of the particle diam- 
eter, the concept of a smooth local interface, and hence 
CWT, breaks down. Naturally, much effort has been de- 
voted to extend CWT, in order to yield a valid descrip- 
tion at short wavelengths&iSaii. The task is challeng- 
ing, because a correct theory must interpolate between 
the continuum (long wavelength) CWT limit, to the mi- 
croscopic (single-particle) limit at short wavelengths. In 




FIG. 1: Schematic representation of an interface within the 
framework of capillary wave theory. The mesh shows the lo- 
cal interface position z — h(x,y), which is assumed to be 
a smooth function of the lateral coordinates x and y. The 
surface h(x, y) then separates phase A (shaded lower region) 
from phase B (upper region). 



simulations and experiments, the situation is slightly less 
problematic. Here, model interfaces are readily prepared, 
and their properties can be probed relatively easily, on a 
wide range of different length scales. 

In this work, Monte Carlo simulations are used to 
study capillary waves in a simple model of a colloid- 
polymer mixture. Our first aim is to test the long wave- 
length predictions of CWT. One of these predictions re- 
lates the capillary spectrum to the interfacial tension. 
The CWT estimate of the interfacial tension, is then com- 
pared to independent estimates, obtained by us in previ- 
ous work. Next, we consider the short wavelength limit of 
the capillary spectrum. In this limit, we observe that the 
capillary spectrum approaches the static structure factor 
of the bulk liquid. Moreover, we find that the transition 
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from the long wavelength CWT limit of the spectrum, to 
the short wavelength limit, is smooth. Our findings are 
then discussed in light of the Mecke-Dietrich theory^ for 
interfaces, and the recent experiment of Ref. 0. 

The outline of this paper is as follows. We first discuss 
CWT. Next, we introduce the colloid-polymer model, 
and describe the simulation method. We then describe 
how the local interface position is extracted from raw sim- 
ulation data. Finally, we present our results, and finish 
with a number of conclusions in the last section. 



II. CAPILLARY WAVE THEORY 

As stated before, the key assumption in CWT is the 
notion of a smooth interface, see Fig. ^ The local in- 
terface position h{x, y) is written as a function of the 
lateral coordinates x and y. By lateral, we mean those 
directions parallel to the plane containing the interface. 
The other important direction is the perpendicular di- 
rection. In this work, the lateral coordinates (x, y) are 
restricted to the interval [0, L], with L the lateral dimen- 
sion of the system. The perpendicular coordinate z is 
restricted to [0, D], with D the perpendicular dimension 
of the system. 



A. Capillary spectrum 



over all pairs of integers (n, m), with n > and m > 0, 
excluding the pair (0,0). The Fourier amplitudes are 
given by 
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h(x,y) sin(x ■ q nm ) dx dy. (5) 



The capillary wave Hamiltonian of Eq. (|2J) now takes the 
quadratic form 
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i.e. an infinite set of decoupled harmonic oscillators, and 
we can use the equipartition theorem to obtain the cor- 
responding expectation values 
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where we introduced C 2 = [A 2 + B 2 )/2, and dropped 
the subscripts nm. Here, q denotes the magnitude of the 
wavevector, (3 = 1/k^T, fee the Boltzmann constant, and 
T the temperature. Note that \Vh\ <C 1 is assumed. One 
should therefore interpret Eq.Q as the limiting form of 
the capillary spectrum for q — > 0. 



The first step in deriving the capillary wave spectrum, 
is to write down the energy cost TL cvl of having a non- 
flat (or undulated) interface with respect to a perfectly 
flat one. Here, the subscript "cw" emphasizes that the 
energy cost is due to capillary waves. The energy can be 
expressed as the product of the interfacial tension 7 and 
the excess area AA. The excess area is defined as the 
area of the undulated surface h(x,y) minus the area of a 
perfectly flat interface, leading to 21 i 22 i 23 i 24 
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with gradient operator V = (d/dx,d/dy). As a first 
approximation, we assume |V/i| <C 1, and thus 
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Vh ■ Vh dx dy, 



(2) 



which defines the capillary wave Hamiltonian. Next, the 
local interface position is expressed as a two-dimensional 
Fourier series 

h(x,y) = A Q0 (3) 

~1~ ^ A nnl COs(x • qnrn) ^nm sin(x • Qmn)' 

n,m 

with x = (x,y) and wavevector q nm = (27m/L, 2iTm/L). 
In the above and following equations, the summation is 



B. Interface broadening 

An important quantity in the investigation of inter- 
faces, is the average density profile p(z), measured in 
a direction perpendicular to the interface. In standard 
mean-field theory, p(z) is given by21*SSi24 



1 \ PA + PB t PA - PB . , 

p{z) = + tanh 



z- z 
Wo 



(8) 



Here, pa (pb) denotes the bulk density of phase A (B), zq 
the location of the interface, and Wq the intrinsic width. 
In this approach, the interface is assumed to be perfectly 
flat, ignoring therefore capillary waves. 

Clearly, in order to describe capillary waves, mean- 
field theory alone is not sufficient. In the presence of 
capillary waves, the interface is no longer flat, but pos- 
sesses a finite mean-squared width. The mean-squared 
width of the local interface position is defined as W 2 W = 
(h 2 ) — (ft,) 2 . If we assume that the modes are decoupled, 
we obtain 
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where Eq. (JTJ) was also used. The mean-squared width 
thus diverges, both in the limit A m j n — > and A max — > oo. 
In a finite system, the magnitude of the smallest wavevec- 
tor is set by the lateral system size: A m i n = 2n/L. Simi- 
larly, one assumes that below some coarse graining length 
a, the notion of a local interface position h(x, y) breaks 
down: A max = 2n/a. The fact that the notion of a local 
interface position on too small scales makes no sense, is 
easily recognized in the example of the d = 2 Ising model 
near its critical pointi Of course, one can clearly define 
the contour separating the region where up spins perco- 
late, from the region where down spins percolate. How- 
ever, the contour contains many overhangs and hence is 
not a single valued function. For the capillary width we 
thus write 

' ML/a), (12) 



W 2 = 

r r cw 
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where one must keep in mind that a is not precisely 
known. 

In a realistic interface at finite temperature, the in- 
trinsic width Wq and the capillary contribution W cw are 
combined. One way to incorporate this into Eq.JHJ, is via 
the convolution approximationi£*2». Here, one assumes 
that the functional form of the density profile in the pres- 
ence of capillary waves is still given by Eq.lJBJ, but with 
a broadened interfacial width W given by 
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In the above equation, the intrinsic width Wo is entan- 
gled with the cut-off contribution a. As a consequence, 
this equation cannot be used to extract Wq from profiles 
obtained in simulations or experiments 15 . In simulations, 
however, Ea. H13|) is still useful because, by varying L, the 
interfacial tension can be measured. 



C. Beyond capillary wave theory 

The discussion of interface broadening is essentially 
based on Eq.JTJ), valid only in the limit q — ► 0. Therefore, 
a cross-over wavevector qg can be identified, above which 
the spectrum is no longer adequately described by Eq.Q. 
In this limit, higher order terms, involving for example 
the bending rigidity of the interface, become important. 
For even larger wavevectors, another cross-over wavevec- 
tor <7t may be identified. It corresponds to wavelengths 
that probe the interface on single particle resolution and 
smaller. Clearly, on such short length scales, the concept 
of a smooth local interface position h(x,y) breaks down. 
We thus identify the following limits: 

1. Long wavelength limit q < q$. The interface is well 
described by a local interface position h(x,y) and 
the capillary spectrum accurately follows Eq.Q- 

2. Medium wavelength limit qs < q < <?t- The inter- 
face is still described by a local interface position, 
but corrections to Eq.Q are important. 



3. Short wavelength limit q > gr- Complete break- 
down of the local interface position concept. The 
microscopic (single-particle) nature of the interface 
cannot be ignored. 

There is consensus regarding the long wavelength 
limit. For example, both Eq.Q and the predicted 
interface broadening, are quantitatively confirmed by 
simulationsi2ii£i2iii2§i2!. Quite the reverse is true in the 
medium wavelength limit. Here, corrections to Eq.Q 
become important, usually quantified in terms of the in- 
dependent interfacial tension 
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The deviations from Eq.Q at large q, now show up as 
deviations of 7(g) from a constant. Note that 7(g) is 
also accessible in experiments, using X-ray or neutron 
scattering at grazing incidence, which probe the height- 
height correlation function^. 

One can conceive 7(g) as an expansion in terms of q. 
According to the Helfrich HamiltonianSS, the next order 
term involves the bending rigidity k, and one obtains 



7(9) =-/o + nq 2 +0(q 4 ), 



(15) 



with 70 the macroscopic interfacial tension. There is 
much controversy regarding the sign of k. For simple 
fluids, Mecke and Dietrich^ developed a theory for 7(g). 
The theory predicts an initial decrease of 7(g), followed 
by a sharp increase of the form nq 2 , implying a positive 
bending rigidity k at large q. As a result, 7(g) contains a 
minimum. Note that for k > 0, the integration in Ea.flOjl 
can formally be extended to infinity, thereby eliminating 
A max from the theory. In recent experiments^^, a min- 
imum in 7(g) was indeed observed, commonly regarded 
as evidence for the Mecke-Dietrich theory. 

The above state of affairs, however, is by no 
means satisfactory. For example, recent Monte Carlo 
simulations 30 ' 31 did not produce any evidence for the ex- 
istence of a minimum in 7(g). Instead, a monotonic decay 
in 7(g) over one order of magnitude is observed. Simu- 
lations thus yield negative values for the bending rigid- 
ity, see also Ref. 0, as would be expected on theoretical 
grounds for systems with short-ranged interactions 3 ^. 

In order to better understand the capillary spectrum, 
it is essential to also consider the short wavelength limit. 
In this limit, the interface is probed on single particle 
length scales and smaller, and will be very rough as a 
result. The concept of a smooth local interface position 
thus breaks down. For theories based on the concept of a 
smooth interface, it will be difficult to describe the tran- 
sition to the short wavelength limit because the micro- 
scopic structure of the fluid becomes important then. In 
experiments and simulations, the situation is less prob- 
lematic: one can simply decrease the wavelength, and 
record how the spectrum changes. 

This approach was tried to some extent by Steckii 3 ., 
who noted that certain features in the capillary spec- 
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TABLE I: Bulk properties of the AO model with R p /R c = 0.8 
for several values of rf p . Listed are the colloid and polymer 
packing fractions in the colloidal vapor and the colloidal liquid 
phase, as well as the interfacial tension 7. All values were 
obtained using histogram reweighting— . 
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trum at large q, are closely related to the structure fac- 
tor S(q) of the bulk liquid. This is to be expected. On 
short length scales, the interface, an essentially macro- 
scopic object, cannot be probed. Instead, one picks up 
the most dominant bulk fluctuations in that case, which 
for a fluid-vapor interface stem from the liquid. The 
simulations presented in this work confirm this picture. 
More precisely, we find that, beyond the first peak in the 
structure factor, the capillary wave amplitudes and the 
structure factor are essentially the same: (C 2 ) ~ S(q). 
Note that S(q) approaches unity, implying 7(5) oc 1/q 2 
in the limit q — > 00. 

Furthermore, our simulations indicate that the bounds, 
qs and <7t, are not sharp but somewhat arbitrary. We find 
that 7(9) evolves smoothly from its constant value 70 at 
low q, to its limiting form 1/q 2 at large q, where (C 2 ) es- 
sentially follows the structure factor. We thus conclude 
that, only in the long wavelength limit, does 7(g) reflect 
a true interfacial property. At short wavelengths, 7(g) is 
largely determined by bulk fluctuations. Since the tran- 
sition is smooth, this implies that in the medium wave- 
length limit, 7(g) is governed by a combination of bulk 
and interfacial fluctuations. 



III. MODEL AND SIMULATION METHOD 

In this work, the phenomena described in the previ- 
ous section are investigated in a model colloid-polymer 
mixture, using grand canonical Monte Carlo simulations. 
Since it recently became possible to visualize the capil- 
lary waves in such systems directlyi&, a simulation in this 
direction seems the most rewarding. Note that colloid- 
polymer interfaces are analogous to fluid-vapor inter- 
faces, albeit on a much larger length scale. A compar- 
ison of our findings, to the abundant literature on the 
fluid- vapor interface in atomic systems, is therefore still 
warranted. 

To describe the particle interactions, we use the so- 
called Asakura-Oosawa (AO) model 3 ^ 3 ^. In this model, 
colloids and polymers are treated as spheres with respec- 
tive radii R c and R p . Hard sphere interactions are as- 
sumed between colloid-colloid and colloid-polymer pairs, 



while polymer-polymer pairs can interpenetrate freely. 
Since all allowed configurations have zero potential en- 
ergy, the temperature plays a trivial role. Throughout 
this work, we consider a size ratio R p /R c = 0.8, and put 
R c = 1 to set the length scale. For R p /R c = 0.8, the AO 
model phase separates into a colloid dense (liquid) and 
colloid poor (vapor) phase, provided the fugacity of the 
polymers z p is sufficiently high. Following convention, we 
use the parameter rf p = z p {Ai:R p /2>) to express the poly- 
mer fugacity, rather than z p itself. To make the analogy 
to atomic systems, rf p may be regarded as the inverse 
temperature. We also introduce the packing fractions 
rji = (4:irRf/3)pi, with pi — Ni/V the number density 
of species i € (c, p), and V the volume of the simulation 
box. 

In Ref. the unmixing behavior of the AO model 
with R v /R c — 0.8 was studied using histogram reweight- 
ing in the grand canonical ensemble. The binodal can be 
found in the same reference, the properties of the state- 
points relevant to this work are listed in Table [I] Inside 
the phase-separated region, two coexisting bulk phases 
separated by an interface can be observed. It is in this 
region where the predictions of CWT can be tested, and 
where the simulations of this work are consequently car- 
ried out. More specifically, we consider state-points with 
an overall packing fraction rj c = 0.134, and rf p ranging 
from 0.9 to 1.2, which is well away from the critical point. 

The simulations of the AO model are performed in a 
rectangular box spanned by the vectors Lx, Ly and Dz, 
using periodic boundary conditions in all three directions. 
Here, x , y and z are unit vectors in the x, y and z direc- 
tion, respectively. In the following, we always use D > i, 
so the interface with the smallest area (and thus lowest 
free energy) is oriented perpendicular to the elongated z 
direction. By using an elongated box, we make sure that 
the interface forms parallel to the (xy)-plane, such that 
L can safely be regarded as the lateral length scale, in ac- 
cord with Fig.^ Note that because of periodic boundary 
conditions two interfaces will actually form. 

To simulate the AO model, we first choose a state- 
point of interest somewhere in the phase-separated re- 
gion of the phase diagram. Next, N c colloidal particles 
are randomly distributed around the center of the simu- 
lation box, with the constraint that colloid-colloid over- 
laps are forbidden. This leads to a bare colloidal system 
without any polymers. The polymers enter the system 
by virtue of a recently developed grand canonical cluster 
mov o 33 i 36 . The characteristic feature of the grand canon- 
ical ensemble is that the number of particles in the system 
is a fluctuating quantity. The idea of the cluster move is 
to insert a number of polymers for each colloid that is re- 
moved, and vice versa. The polymers thus enter the simu- 
lation box via a repeated application of cluster moves. To 
keep the colloid packing fraction fixed at rj c — 0.134 (and 
hence stabilize the interface) , we additionally restrict the 
cluster moves such that only configurations with N c and 
N c + 1 colloids are accepted. No constraint is put on the 
number of polymers though. Note that the acceptance 
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FIG. 2: Schematic representation of a phase-separated 
colloid-polymer mixture. Colloidal particles are drawn as cir- 
cles, the polymers are not shown. The colloid dense phase 
is centered around the middle of the box at z = D/2. The 
shaded regions represent interfaces, located around Z\ and z%, 
and respective widths wi and u>2- See also details in text. 



probability of the cluster move depends on the fugacity of 
both the colloids and the polymers^, and this is where 
1]^ comes into play. The colloid fugacity is tuned such 
that the states with N c and N c + 1 are visited equally 
often on average. 

This procedure provides a rather unbiased way to 
study phase separation. By starting with a random con- 
figuration of colloids, we make sure that any observed 
phase-separation is thermodynamically driven, and not 
some simulation artifact. Moreover, we do not specify 
the coexistence densities of the phases in any way. In- 
stead, these densities are an output of the simulation, and 
should agree with the densities of the bulk phase diagram 
obtained, for instance, from histogram reweighting. This 
is used as a test for equilibration. After equilibration, 
we continue to simulate until around 100 uncorrelated 
configurations have been obtained. The resulting con- 
figurations are then used for the interface analysis to be 
discussed next. During the simulations, the grand canon- 
ical cluster moves are used in conjunction with random 
displacements of single particles, typically attempted in 
a ratio (5:1), respectively. 



IV. INTERFACE EXTRACTION 

By using the method outlined in the previous sec- 
tion, large numbers of phase-separated colloid-polymer 
configurations can be generated. In this section, we 
explain how the local interface position h(x,y) is ex- 
tracted from these configurations, using the block anal- 
ysis methodise. This method was used successfully be- 
fore in the study of interfaces in polymer blends. 
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FIG. 3: Profile of the colloid packing fraction (solid line) and 
polymer packing fraction (dashed line) along the elongated z 
direction at -q p = 1.1. The profiles were obtained from a sin- 
gle configuration of size {L — 60; D = 120}. This corresponds 
to approximately 13,000 colloidal particles and 130,000 poly- 
mers. The number of slabs used in calculating the profiles 
was N = 100. The horizontal dotted lines represent the bulk 
packing fractions obtained from histogram reweighting taken 
from TableQ] Arrows mark the bounds of the interface regions 
obtained by fitting to Eq.JSJ. Each interface region contains 
around 2000 colloids and 16,000 polymers. 



A. Global interface localization 

Before applying the block analysis method, the inter- 
faces in the configuration are isolated from the bulk re- 
gions. In doing so, we ensure that the local interface posi- 
tion h(x, y) is not influenced by density fluctuations that 
occur deep inside the bulk. A schematic representation of 
a phase-separated configuration is shown in Fig. |5J The 
colloids are drawn as circles, the polymers are for clarity 
not shown. 

The use of periodic boundary conditions leads to the 
formation of two interfaces. The interface regions are 
shaded gray in Fig. One interface is located around z\ , 
the other around z%. The respective widths are w\ and 
W2 ■ To obtain numerical values for the interface locations 
and the widths, we divide the box into N equal slabs 
perpendicular to the z direction. Each slab has area L 2 
and width Az = D/N. The colloid density profile along 
the z direction is now estimated by p{z) — nj (AzL 2 ), 
with n the number of colloids in the slab centered around 
z. A similar expression holds for the polymer density 
profile. 

Typical profiles for the AO model obtained in this way 
are shown in Fig. |3 with density converted to packing 
fraction. The solid line shows the colloid packing fraction 
as a function of z, the dashed line shows the correspond- 
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FIG. 4: Schematic illustration of the block analysis method. 
The interface region is divided into rectangular segments. In 
each segment, the Gibbs surface defines a local interface po- 
sition zq. See details in text. 



FIG. 5: Schematic density profile (colloid or polymer) along 
the z direction in a single segment. The Gibbs surface is 
located at zg, such that areas A and B are equal. 



ing polymer packing fraction. Important in these profiles 
is the height of the plateaus, which must be in agreement 
with the bulk packing fractions obtained from histogram 
reweighting. The latter values are listed in Table [IJ and 
represented in Fig. CUby the horizontal lines. We observe 
that the plateaus are in good agreement with Tabled in- 
dicating that the configuration has reached thermal equi- 
librium. 

Next, we fit the hyperbolic tangent of Eq.JSJ to the 
profiles of Fig. [3] We emphasize that these are two- 
parameter fits in the width and interface location only. 
For the coexistence densities, one uses the values of the 
bulk phase diagram. To measure the location z\ of the 
left interface and its width w\ , the range of the fit is set 
to < z < D/2. For the right interface, the appropriate 
range is D/2 < z < D. One still has the choice to fit 
to the colloid profile, or to the polymer profile. Since 
the number of polymers is much larger than the number 
of colloids, a fit to the polymer profile seems to be the 
most sensible. To ensure that the entire interface region 
is captured, we multiply the width W obtained from the 
fit by a factor a > 1. The bounds of the interface region 
are thus written as: 



z± = z ± aW. (16) 



As an illustration, the best fit parameters for the poly- 
mer profile in Fig. are {z = 36.5; W = 2.37} (left) 
and {zq = 83.5; W = 2.86} (right). The corresponding 
bounds are marked by the arrows in Fig. [31 where we 
used a = 3. 



B. Block analysis 

Having identified the interface regions in a single con- 
figuration, the next task is to extract the local inter- 
face position h(x,y). For this we use the block analy- 
sis method, which is schematically illustrated in Fig. 01 
The figure shows one of the interface regions, and con- 
tains all those particles whose z coordinates are between 
z_ < z < z + . Note that in Fig. 01 the elongated z 
direction points into the plane of the paper. The inter- 
face region is split up into rectangular segments of size 
(L/ub) x (L/ub) x {z+ — z~). Here, the parameter «b 
is an integer called the block factor. The total number 
of segments thus equals n|. An example of one such 
segment is also shown in Fig. 01 

We now assign a local interface position zq in each 
segment, in spirit of the Gibbs surface. To this end, we 
assume that the colloid (polymer) density profile p s (z), 
along the z direction in a single segment, is similar to 
Eq.lJSJ. The Gibbs surface runs perpendicular to the z 
direction, such that the shaded regions in Fig. have 
equal area. For the areas A and B we may write 

/•ZG 

A= p s {z)&z - pa{zg - z-), (17) 

J Z- 

f z+ 

B = p B {z+ - zg) - I Ps(z)dz, (18) 

J ZG 

with pa and ps the density of the colloids (polymers) 
in the vapor and liquid phase, respectively. Equating 
the above two expressions, and multiplying by the lateral 
area of a single segment, we obtain 

n(n B /L) 2 = pa(z g - z_) + p B {z+ - z G ), (19) 
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FIG. 7: Capillary amplitudes, calculated using Eq.@ and 
Eq.@, for several values of the block factor tib. The data 
were obtained in simulations of the AO model at if p = 1.2 and 
box size {L — 60; D — 120}. Averages over approximately 
100 configurations were taken. The lines at low q represent 
the CWT prediction of Eq.Q, with 7 taken from Tableland 
no other free parameters. For clarity, each data set has been 
multiplied by tib- 



FIG. 6: Influence of the block factor na on the local interface 
position h(x, y). Each of the above h(x, y) was extracted from 
the same AO configuration at rf p = 1.1 and box size {L = 
60; D = 120}. The top frame shows h(x,y) obtained using 
tib — 8, and the lower frame using ub = 24. 



with n the number of colloids (polymers) in the segment. 
This equation is easily solved for zq, because pa and ps 
are known from the phase diagram. One attractive fea- 
ture of Ea. (|19fl is that the precise form of the density 
profile need not be specified. This is also required, if 
the short wavelength limit is to be probed. Large block 
factors iib must then be used, and the number of par- 
ticles in a single segment will predominantly be zero or 
one. Obviously, with such low numbers, a smooth pro- 
file cannot be defined in any case. We emphasize that 
by using Ea. ((T§l) . two estimates for zq may be obtained: 
one determined by the number of colloids, and one by 
the number of polymers. Note also that zq may well be 
located outside the region bounded by z_ and z + . 



C. Influence of parameters 

The local interface position, as defined in the previous 
section, thus depends on three adjustable parameters, 
namely: (1) the block factor Ub, (2) a in Ea. (|16fl . and (3) 
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FIG. 8: Capillary amplitudes obtained using several values 
of a in Ea. (ll61 as indicated. The straight line represents the 
CWT form of Eq.0, with 7 taken from TableQ The spectra 
were measured at rf v — 1.2 and box size {L = 60; D = 120}. 
For clarity, the data have been smoothed. 
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FIG. 9: Capillary amplitudes obtained using the polymers 
(triangles) and the colloids (circles) in Eg. 1191 . The straight 
line is the CWT form of Eq.Q, with 7 taken from Table □ 
The spectra were measured at rf p = 1.2 and box size {L = 
60;L» = 120}. 



FIG. 10: Correlation time r of the amplitudes A„o, as a 
function of the wave number n, for several system sizes as 
indicated. The correlation time is expressed in seconds of 
CPU time on a 2.4 GHz Pentium machine. The local interface 
position h(x, y) was extracted from AO configurations with 
rip — 1.1 and using block factor «b = 20. 



whether n in Eq. 119fl is taken to be the number of colloids 
or the number of polymers. It is important to determine 
the sensitivity of the results to these parameters. 

The influence of the block factor is discussed first. 
Here, we use a = 3, and n in Ea. l|19|) is taken to be 
the number of colloids. The local interface according 
to our method is not defined in continuous space, but 
on a grid of lattice spacing L/ub, see Fig. 0J The in- 
dices (n,m) in the summation of Eq.@ can therefore 
run up to tib/2 only. As a result, the shortest wavelength 
that can be sampled equals a — 2L/ub 1 or alternatively 
A max = ttub/L. One may therefore conceive Hb as a fil- 
ter: via ne, all capillary modes with q > A max are effec- 
tively filtered out. The effect in real-space is illustrated 
in the snapshots of Fig. The AO configuration from 
which the local interface position is extracted, is the same 
in both snapshots, but ub is varied. For small ub, only 
long-wavelength variations are visible. For larger values, 
short-wavelength features become visible, and the inter- 
face appears rougher. The effect of tib in Fourier-space 
is shown in Fig. [7| Here, we measured the Fourier am- 
plitudes as a function of q, for several values of tib- By 
increasing iiBi the range of the data can be extended to 
larger values of q, but otherwise the spectra are unaf- 
fected. In particular, note that the spectra for q — > 
are well described by the CWT prediction of Eq. J7J) , ir- 
respective Of 7T-B- 

Next, a is varied at fixed block factor tib — 64, with n 
in Eq. l|19|l again being the number of colloids. This result 
is shown in Fig. |SJ For a = 5, the bounds z± lie deep 
inside the bulk, so considering larger values is not sensi- 



ble. For small q, we observe that the spectra are rather 
insensitive to a and, moreover, in good agreement with 
Eq.Q. For q sa 0.2 and above, systematic deviations 
become visible. 

Finally, Fig. shows the effect of using the polymers 
to determine zq in Eq . (|19|) . instead of the colloids. Here, 
ub = 64 and a — 3 were used. The agreement with 
Eq.Q is again demonstrated for q — > 0, but at higher 
values deviations become visible. Note that L 2 (C 2 ) de- 
viates from the CWT behavior at low q towards larger 
values. This already indicates that 7(g) < 7(0) for all 
nonzero q, see Eg. (11411 . 

In summary, the block factor tib does not influence the 
capillary spectrum in any significant way; it merely sets 
the range in q that can be sampled. In contrast, a in 
Ea. (|16|) . and n in Ea.Q19[l. do affect the spectrum, but 
only at relatively large values of q. 



V. RESULTS 

We now proceed with presenting our main results, 
whereby the long and short wavelength limit of the cap- 
illary wave spectrum are treated separately. 



A. Long wavelength limit: g< 1 

As demonstrated before, the capillary spectrum in the 
long wavelength limit, does not depend sensitively on the 
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FIG. 11: Capillary spectrum of the AO model in the long 
wavelength limit, for several values of rf p as indicated (note 
the double logarithmic scale). Lines represent the CWT form 
given by Eq.©. All spectra were obtained using box size 
{L = 60; D = 120} and block factor n B = 100. 



FIG. 12: Capillary width as a function of the block fac- 
tor (note the logarithmic jib scale). The data were obtained 
using box size {L = 60; D = 120}. Averages over approx- 
imately 100 configurations were taken. The lines represent 
Eg. 1201 . with 7 taken from TableQ and C obtained by fitting 
to the lowest few ns values. 



parameters used to extract the local interface position. In 
this section, we use a = 3 in Eq. I|16l) , and the local inter- 
face position zg in Ea. (|19|) . is determined via the colloids. 
To properly model long wavelength capillary modes, large 
system sizes are required. An important issue is then the 
correlation time r of these modes. It is well-known that 
r increases strongly with the size of the simulation box, 
such that extensive simulations are required^. For the 
AO model, this is illustrated in Fig.^J| where r is plot- 
ted, for several system sizes. The data were obtained 
by measuring the amplitudes A n o(t), given by Eq.Q, 
at evenly spaced times t, with wave number n > 1. 
From the decay of the autocorrelation function, r was 
obtaine d 37 ! 38 . As expected, r decreases with the wave 
number. In a simulation, the important quantity is the 
correlation time of the amplitude with the smallest wave 
number: A\q. It is crucial to ensure that the duration of 
the simulation exceeds the time that governs the decay 
of Aiq. Fig. 1101 shows that the corresponding correla- 
tion time increases strongly with the size of the system. 
We found that a system size of {L = 60; D = 120} is 
sufficient to perform a long wavelength analysis of the 
interface. Since around 100 uncorrelated configurations 
are required to obtain the spectrum with reasonable ac- 
curacy, this amounts to approximately four days of CPU 
time. 

The capillary amplitudes are shown in Fig. 1111 for 
several values of rfL. Also included in the figure is the 
CWT form of Eq.0, with 7 taken from Table [fl and 
no other adjustable parameters. We observe that CWT 
is confirmed for all values of rf p considered by us, up to 
qs ~ 0.2, corresponding to a length scale of around 15 



colloid diameters. 

Next, interface broadening is considered. First, note 
that the logarithmic dependence of W on L, see Ea. lllHI) . 
was already confirmed by us in Ref.[l2 (for rf p — 1.1), but 
with a factor of two discrepancy in the definition of the 
width. This yields an incorrect estimate of the interfacial 
tension. If, however, this error is accounted for, a fit to 
the data yields 7 = 0.109 ± 0.004. The deviation from 
histogram reweighting is around four percent, see TableQ] 

An alternative method to study interface broadening, 
is to vary the coarse graining length a at fixed L, and 
use Eq. (|12|l to obtain the interfacial tension. To use this 
method, a single set of configurations, obtained for one 
value of L, is sufficient. For the coarse graining length, 
we may write a oc 1/ub, so Eq. (|12I) becomes: 

W? w = ^-\n(n B /C) 7 (20) 

with 71b the block factor and C a constant. In this ap- 
proach, the local interface position h(x, y) needs to be ex- 
plicitly calculated. Per interface region, n| samples hi of 
the local interface position are obtained. The variance in 
these samples yields an estimate for W C 2 W = (/if) — (hi) 2 , 
which is then averaged over the different configurations. 
In Fig. El w e show the width W 2 W , as a function of ub, 
for a number of different rf p . The figure strikingly il- 
lustrates the broadening of the interface as the grid, on 
which the interface is studied, becomes finer, and more 
capillary modes are taken into account. For small np, the 
broadening is in good agreement with Eq. H20|) , shown by 
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FIG. 13: Capillary amplitudes for several values of a as in- 
dicated. Here, the colloids were used in Eg. 11911 to deter- 
mine zg- The line at low q is the CWT prediction of Eq.JTJ. 
The data were obtained in simulations of the AO model at 
?7p = 1.0, using box size {L = 20, D — 40}. Also shown is the 
static structure factor S(q) of the bulk colloidal liquid. Note 
the double logarithmic scale. 



the lines in Fig. 1121 Here, 7 was taken from Tabled and 
C was obtained by fitting. For large ne, deviations be- 
come visible. The range of validity varies from ne ~ 4 
(rf = 0.9) to riB ~ 8 (r/p = 1.2). The corresponding 
wavelength is approximately 8—15 colloid diameters, or 
equivalently q$ » 0.2 — 0.3. 



B. Short wavelength limit: g> 1 

We now turn to the short wavelength limit of the cap- 
illary spectrum. In this case, both the parameter a in 
Eci. (|16fl . as well as the precise choice (colloids or poly- 
mers) for n in Eq. l|19|l . become important. To illustrate 
the effect, we show in Fig. 1131 and Fig. ^1 the capillary 
amplitudes for several values of a, at 77^ = 1.0. The data 
in Fig. 1131 were obtained using the Gibbs surface defined 
by the colloids as local interface position. In Fig.^J the 
Gibbs surface defined by the polymers was used. In both 
figures, the line at small q represents the CWT form of 
Eq.0. Also shown in Fig. ^3 is the static structure 
factor S(q) of the bulk colloidal liquid (multiplied by a 
factor of 30). 

One important observation is that, for a sufficiently 
high, the local Gibbs surface defined by the colloids and 
the polymers, both reproduce the expected CWT result 
for q — ► 0. Note also that, compared to Fig. ^2 the 
system size is now much smaller, so the agreement with 
Eq.Q is not as convincingly demonstrated as before. 

More importantly, we observe in Fig. 1131 that for 



FIG. 14: The analogue of Fig. 1131 but using the polymers in 
Eg. 1191 to extract the local interface position. 



gx ~ 1.2 and above, the capillary wave amplitudes es- 
sentially follow the static structure factor. This illus- 
trates the point made earlier, namely that, for high q, the 
capillary spectrum is determined by the most dominant 
bulk fluctuations. The same effect is visible in Fig. 1141 
where we recall that, in the AO model, the polymers are 
ideal (hence no structure at large q). Furthermore, we 
deduce from Fig. 1131 and Fig.^J that the limiting value 
lim^oo L 2 {C 2 ) is proportional to a: in other words, pro- 
portional to the width of the slabs in Fig. [2J or conse- 
quently, the number of particles involved in calculating 
L 2 (C 2 ). This is again consistent with the view that the 
capillary amplitudes approach the (unnormalized) static 
structure factor of the bulk (note that the static struc- 
ture factor is usually normalized by dividing through the 
total number of particles). 

Finally, we observe that the capillary amplitudes 
evolve smoothly from L 2 (C 2 ) cx 1/q 2 at low q, 
to L 2 (C 2 ) cx S(q) at high q. This implies, for 
the q-dependent interfacial tension, lim g _ (00 7(g) cx 
1/q 2 , see Fig. In agreement with (most) other 

simulationsi^SiSi, we observe a strong reduction in j(q) 
with increasing q, but no sign of a minimum. 



VI. DISCUSSION AND CONCLUSIONS 

First, we note that our data in the long wavelength 
limit, are consistent with CWT. The interfacial tension, 
obtained from both the spectrum and the broadening of 
the interface, is in good agreement with previous inde- 
pendent estimates^. This result is encouraging, because 
it is difficult to simulate the long wavelength limit accu- 
rately, due to the large system sizes that are required. 
For the AO model, we observe that CWT remains valid 
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FIG. 15: Momentum dependent interfacial tension 7(g) for 
the AO model, as a function of q 2 , for several values of rf p 
as indicated. The data were obtained using box size {L — 
60, D = 120}, with n in Eg. 11911 being the number of colloids, 
and a = 3. 



down to wavelengths of around 10 colloid diameters, in 
reasonable agreement with Ref. The data in Fig. ^] 
and Fig. also show that the deviations become more 
pronounced closer to the critical point (so for low values 
of r/p). This is in agreement with Ref. 39. 

At shorter wavelengths, CWT gradually breaks down. 
We observe a sharp reduction in 7(9), followed by an 
asymptotic decay of the form 7(g) oc l/q 2 . In the asymp- 
totic regime, 7(g) is essentially determined by the static 
structure factor, see Fig. ED and Fig. ^] The form of 
7(g) shown in Fig. 1151 is consistent with most other 
simulationsi^22i2i (an exception is Ref. |3^, in which an 
entirely different definition of the local interface position 
is used). 

However, the problem involved in our analysis (as well 
as in related previous work) is that the concept of a local 
interface in terms of the Gibbs surface becomes ques- 
tionable as the considered lateral length scale becomes 
smaller and smaller. Since in the theory of Mecke and 
Dietrich^ a somewhat different approach is used, it is 
not clear if there is a real contradiction between our re- 



sults and their treatment (note also that we work with a 
model exhibiting strictly short range forces, while Ref. 
employs long range van der Waals forces). Our results 
clearly point towards the need of a different analysis of 
simulations (and experiments) suitable to consider the 
interplay of bulk and interfacial fluctuations on inter- 
mediate and short length scales. Note that the pro- 
nounced size-dependence of the simulated interfacial pro- 
files makes our analysis in terms of the intrinsic interfacial 
width rather meaningless as well, and this complicates a 
direct comparison of our data to theoretical studies^2*^i. 

The problem that the concept of a local interface be- 
comes ill-defined on the microscopic scale (correlation 
lenght in the bulk, or interparticle distances in the fluid 
in the extreme case) also calls into question the extent to 
which the concept of a wavevector-dependent interfacial 
tension 7(g) is really meaningful. Fig. ^] and Fig. ITU 
indicate that one must be careful with the choice of the 
parameter a, and only for a > 1 do we observe con- 
vergence of the data for small q. At large g, however, 
the dependence on a is very pronounced (reflecting es- 
sentially the number of particles used in the calculation) . 
Nevertheless, whatever the value of a, the data always 
follow the static structure factor in the bulk, see Fig.lT3l 
This finding is sensible because if one probes the inter- 
face on such small scales, the local structure reflecting 
how the particles in the fluid are stacked, should show 
up in the analysis. 

It remains a challenge for future work to try to extract 
from simulations the precise analog of the structure factor 
obtained in scattering experiments^*^ and to conduct 
such an analysis with the present data. 
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